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A  characteristic  mmdimensional  distortion  number  N^i  derived  in  the  1970s  that 
allows  inference  of  the  degree  of  nonlinear  refraction  or  thermai  blooming  associated  with 
an  atmospheric  laser  path.  For  a  continuoi4S-wave  ( CW)  laser  with  a  Gaussian-shaped 
beam,  the  distortion  number  is  a  function  of  several  variables  including  laser  power  and 
aperture  size,  optical  wavelength,  atmospheric  absorption  and  extinction,  index  of 
refraction,  temperature,  air  density,  and  the  air  speed  or  flow  transverse  to  the  laser  beam. 
Scenario-dependent  calculations  of  atmospheric  distortion  number  N^f  are  developed  for 
different  geographic  regions  and  seasons  using  the  Air  Force  Research  Laboratory's 
global  thermosonde  database,  the  HITRAN  molecular  spectroscopic  database,  and  global 
climatological  aerosol  model  extinction  profiles.  Tactical  air-to-ground  scenarios  are 
described  as  a  function  of  altitude,  target  distance,  and  taser-to*target  azimuth  angle  for 
the  COIL  wavelength  (IJ15  pm).  The  results  are  interpreted  in  light  of  seasonal  and 
geographical  factors  as  well  as  path-integrated  moisture. 

KEYWORDS:  Disionioti  number*  Molecular  absorpiion*  Rcfraciion-  Thermal  blooming 


Nomenclature 

a  laser  apenure  radius,  m 

turbulent  temperature  structure  constant, 
c>  speciHc  heat  of  air  at  constant  pressure,  J  kg^’  K'* 

4  /o  area-averaged  intensities  at  focal  plane  center  for  blooming,  diffraction-limited 
propagation,  respectively,  W  m”^ 

!p  normalized  peak  inradiance 

k  laser  wave  number, 

absorption  number 
distortion  number 
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Nf  Fresnel  number 

N„.  slew  number 

n  index  of  refraction  of  air  at  laser  wavelength 

F  laser  power,  J  s"* 

laser  beam  radius,  m 

S  atmospheric  path  range  or  total  distance,  m 

SPW  slant  path  precipi  table  water  cm 
j  distance  along  laser  path  from  source  to  receiver  m 

T  virtual  air  temperature.  K 

w  water  vapor  mixing  ratio,  g  kg“' 

Co  nonlinear  interaction  parameter  m 

af  atmospheric  total  absorption  coefficient,  m“' 

€  index  of  refraction  perturbation 

A  laser  wavelength, 

Vff  beam  transverse  air  speed,  m  s' ' 

p  air  density,  kg 

Pti  density  of  dry  air,  kg 

Pf  de nsi ty  o f  1  iq u i d  wate r  kg  m  ' ^ 

<j  atmospheric  extinction,  ni^^ 


1.  Introduction 

Nonlinear  refraction,  often  referred  to  as  “thermal  blooming^  is  the  distortion  of  the  laser 
irradiance  held  caused  by  its  interaction  with  the  mean  index  of  refraction  field  of  the  atmo¬ 
sphere.  The  effects  of  thermal  blooming  are  related  to  the  local  heating  of  the  atmosphere 
caused  by  the  absorption  of  the  high-energy  laser  (HEL)  energy  and  the  advection  of  heat  by 
beam-transverse  flow  a,s  the  laser  wave  front  propagates  thniugh  the  atmosphere.  This  inter¬ 
action  between  the  laser  power  and  the  atmosphere  is  nonlinear.  Thermal  blooming  effects 
can  be  numerically  computed  by  integrating  the  paraxial  wave  equation  approximation  to 
Maxwell's  equation  (e.g..  Refs.  7  and  1 7).  This  equation  includes  a  simplihed  steady-state 
relation  between  refractive  and  irradiance  fields  that  is  deri  ved  from  the  linearized  equations 
of  hydrodynamics, 

Wallace  and  Camac'^  derived  a  nonlinear  interaction  parameter  or  distance  zo*  They 
described  as  the  propagation  distance  at  which  the  nonlinear  refraction  effects  caused  by 
transverse  density  gradients  distributed  over  the  laser  beam  width  will  distort  the  irradiance 
distribution  significantly.  For  a  continuous- wave  (CW)  laser,  Z{\  is  expre,ssed  as  the  ratio  of 
the  laser  beam  radius  to  the  square  root  of  the  normalized  index  of  refraction  perturbation 
e: 

(n  -  \)ct,F 

.  lU 

where  n  is  the  refractive  index  of  air  at  the  laser  wavelength,  or,  is  the  atmospheric  total 
absorption  coefficient  at  the  laser  wavelength,  P  is  laser  power,  p  is  air  density,  Cp  is  the 
specific  heal  of  air  at  constant  pressure,  T  is  the  virtual  air  temperature,  and  is  the  beam 
normal  or  transverse  air  speed.  The  virtual  air  temperature,  which  satisfies  the  equation  of 
state  for  dry  air  in  a  moist  atmosphere,  is  a  function  of  air  temperature  and  water  vapor 
mixing  ratio.  The  form  of  Eq.  ( I )  implies  a  homogeneous  atmosphere.  The  transverse  air 
.speed  is  assumed  to  be  subsonic. 
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Bradley  and  Herrmann’  delined  a  similar  expression,  calling  it  the  nondimensional  dis¬ 
tortion  number  and  used  it  to  quantify  the  degree  of  thermal  blooming  associated  with 
an  atmospheric  path  range  or  focal  length  5.  Gebbardr"^  generalized  their  expression  for 
a  nonhomogeneous  path  [see  Eqs.  (4.22)  and  (4.24)  in  Ref.  5].  For  a  CW  laser,  this  is 
expressed  as 

Nj  =  — ■ —  /  — - - —  djT,  (2) 

a  J<i  Cpn(s)p{s}v„{s) 

where  the  variable  symbols  are  the  same  as  described  in  Eq.  ( I ).  Here  P,  is  the  laser  power 
at  the  laser  aperture  of  radius  a,  ^  is  the  laser  wave  number,  a  is  the  atmospheric  extinction, 
and  s  is  the  distance  along  the  laser  path  from  source  to  receiver 


2.  Relation  of  to  Laser  Beam  Distortion 

The  relation  of  N^i  to  laser  beam  irradiance  field  distonion  due  to  nonlinear  refraction 
is  illustrated  by  integrating  the  paraxial  wave  equation  of  Bradley  and  Herrmann  [Ref.  2. 
Eq.  (3)1  in  the  propagation  direction  for  a  CW  laser  with  circular  aperture.  This  nondimen¬ 
sional  differential  equation  includes  both  the  atmospheric  distonion  number  and  a  set  of 
three  other  nondimensional  numbers  including  the  Fresnel  number  Nf,  absorption  number 
Nfi,  and  slew  number  We  integrated  this  equation  using  the  altemate-direciion  implicit 
finite  element  method  (Ref.  4.  Chapter  4).  Figure  la  depicts  an  image  of  normalized  irradi¬ 
ance  on  the  focal  plane  for  diffraction- limited  effects,  i.e..  accounting  for  no  atmosphere. 
Figure  lb  shows  the  irradiance  field  for  steady-state  uncompensated  propagation  of  the  CW 


a)  DifTraction  limited  (no  atmosphere)  b)  Uncompensated  ihcrmal 

blooming 

Fig.  L  Simulated  images  of  the  focal  plane  irradiance  depicting  a)  diffraction  limited  and 
b)  uncompensated  blooming  occurring  with  steady-state  CW  HEL  Gaussian  beam  prop¬ 
agation  in  a  homogeneous  atmosphere.  The  propagation  direction  is  normal  to  the  page, 
and  the  beam-transverse  flow  is  from  left  to  right  in  panel  b.  The  transverse  coordinates  are 
scaled  by  beam  radius.  Irradiance  levels  are  proportional  to  shade  lightness. 
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Fig.  2.  Normalized  peak  intensity  as  a  function  of  atmospheric  distortion  number 
Even  with  phase  compensation,  significant  intensity  reduction  begins  for  >  ICK).  The 
black  circle  corresponds  to  for  the  uncompensated  blooming  depicted  in  Fig.  lb.  The 
slew  number  No,  and  absorption  number  N^t  are  also  shown  in  the  diagram  of  Bradley  and 
Herrmanns 


laser  through  a  homogeneous  atmosphere  where  the  distortion  number  A/j  exceeds  UK). 
Effects  of  refractive  turbulence  on  the  laser  irradiance  field  are  not  included.  The  transverse 
flow  speed  v„  is  directed  from  left  to  right  in  the  figure.  In  the  presence  of  an  atmosphere 
and  thermal  blooming  (Fig.  lb),  significant  phase  and  amplitude  distortions  of  the  laser 
irradiance  occur. 

Figure  2  depicts  the  change  of  the  normalized  peak  irradiance  ip  as  a  function  of  distonion 
number  for  uncompensated  and  phase- compensated  atmospheric  effects.  Here.  Ip  is 
defined  as 


tp  =  (3) 

h 

where  /*  and  /a  are  the  area-averaged  intensities  near  the  center  of  the  focal  plane  for  thermal 
blooming  and  diffraction-limited  propagation,  respectively.  The  ratio  4//(j  is  similar  to  the 
Strehl  intensity  defined  in  Bom  and  Wolf  IRef.  L  Chapter  9,  Eq,  (9)].  Increasing  while 
holding  the  atmospheric  parameters  fixed  is  equivalent  to  increasing  the  laser  aperture  power. 
In  the  case  of  uncorrected  blooming,  the  normalized  peak  irradiance  Ip  begins  to  decrease 
at  relatively  small  values  of  yVj'^30.  Phase -com  pen  sated  propagation,  using  Eq.  (4)  of 
Bradley  and  Herrmanns  permits  the  normalized  peak  irradiance  to  increase  for  values  of 
Nit  5  ItK),  after  which  Ip  drops  rapidly  once  the  distortion  number  exceeds  a  critical  value. 
Thus,  increasing  the  distortion  number  whether  through  a  power  increase  or  atmospheric 
absorption  increase,  can  eventually  lead  to  decreased  intensity  in  the  focal  plane,  even  when 
the  beam  is  phase  compensated  near  the  aperture. 
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3,  Estimating  the  Atmospheric  Distortion  Number 

The  daia  source  for  the  atmospheric  distortion  number  analyses  are  thermosonde  vertical 
profile  measurements.  Thermosondes  are  balloon-borne  payloads  that  measure  in  situ  the 
turbulent  temperature  structure  constant  Cf  approximately  every  7-8  m  in  the  vertical  using 
fine-wire  probes,  Cf  measurements  and  the  associated  refractive  turbulence  are  not  included 
in  the  distortion  number  calculations  here.  In  addition  to  Cj,  the  thermosondes  measure 
pressure,  temperature,  humidity,  and  horizontal  wind  velocity  using  attached  radiosondes. 
In  practice,  synoptic  radiosonde  data  could  be  substituted  for  thermosonde  data  in  these 
estimates,  but  this  would  result  in  significantly  lower  vertical  resolution. 

Equation  (2)  is  evaluated  by  integrating  along  an  air-to-ground  slant  path  (nonrefracted 
ray)  assuming  a  spherical  earth.  For  the  slant  path,  the  upper  end  is  defined  in  terms  of 
altitude  above  mean  sea  level  (MSL),  the  lower  or  ground  end  is  set  at  5  m  above  the 
surface,  and  the  path  length  is  defined  in  terms  of  earth  surface  distance.  The  calculations 
are  performed  for  a  range  of  altitudes — up  to  6  km  above  MSL — and  for  earth  surface 
distances — 5-35  km — which  bound  approximate  tactical  scenarios.  These  altitudes  and 
surface  di.stances  approximate  a  range  of  tactical  scenarios. 

The  evaluation  of  the  beam  transverse  wind  speed  v„  along  the  path  includes  both  the 
beam  slew  velocity  and  horizontal  wind  velocity  measured  by  the  thermosonde.  The  slew 
speed  path  component  is  evaluated  using  an  upper-end  air  speed  of  100  m/s  with  360-deg 
heading,  a  fixed  (stationary)  point  at  the  ground,  and  three  different  target  azimuths  (0,  45, 
and  90  deg)  relative  to  the  upper-end  heading.  The  laser  aperture  radius  ts  set  at  50  cm, 
and  the  laser  power  is  10^  W  (Ref.  5,  Table  4.2),  The  wavelength  of  1.315  u,sed  for 
the  Chemical  Oxygen- Iodine  Laser  (COIL)'^  h  selected.  Table  I  lists  the  fixed  parameters 
used  in  the  evaluation  of  Eq.  (2).  The  air  density  p,  refractive  index  n,  and  temperature  T 
along  the  path  are  determined  from  the  thermosonde  measurements. 

The  atmospheric  absorption  of  and  extinction  a  at  the  laser  wavelength  are  important 
factors  in  determining  the  distortion  number  for  tactical  scenarios.  However,  their  deter¬ 
mination  Is  not  as  straightforward  as  for  the  other  variables.  In  this  study,  the  absorption 
coefficient  represents  the  sum  of  molecular  line-by-line  absorption,  continuum  absorp¬ 
tion,  and  aerosol  absorption,  which  should  be  treated  differently  from  molecular  absorp¬ 
tion  for  thermal  blooming,  as  described  below.  The  extinction  coefficient  i,s  the  sum  of 


Table  1,  Fixed  parameters  for  calculation  of  the  atmospheric 
distortion  number  Nti  using  Eq,  (2) 


Fixed  parameter  type 

Magnitude 

Laser  beam  diameter 

1  m 

Laser  wavelength 

1.315  ^m 

Laser  power 

lO^  W 

Upper-end  air  speed 

100  m/s 

Lower-end  ground  speed 

Om/s 

Lower-end  altitude 

5  m  AGL 

Upper-end  altitudes 

0.5-6  km  MSL 

Upper-end  heading 

360  deg 

Lower-end  azimuth  angles 

0,  45,  90  deg 
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Fig-  3-  Average  surtace  visibilities  (meteorological  range)  at  differem  locaiions  and  during 
different  seasons  trir  thermosonde  measurement  campaigns  based  on  the  GADS  aer(>so) 
data  sets- 


total  absorption-  molecular  scattering,  and  aerosol  scattering-  The  molecular  absorption 
at  1.315  fim  is  determined  from  line-by-line  calculations  using  the  HITRAN  molecular 
database.'-^  including  the  molecular  continua,  and  the  thermosonde- measured  vertical  pro¬ 
files  of  temperature-  pressure-  and  relative  humidity.  The  aerosol  extinction  and  absorption 
parameters  are  derived  using  thermosonde  data  and  the  Optical  Properties  of  Aerosols  and 
Clouds  (OPAC)/Global  Aerosol  Data  Sets  (GADS).®’^^  GADS  provides  a  climatological 
database  of  aerosol  absorption  and  extinction  profiles  for  a  5  x  5  deg  global  grid  given 
geographical  location-  season  (summer  or  winter)*  and  surface  relative  humidity.  One  of 
seven  possible  aerosol  altitude  profile  types  is  assigned  to  each  grid  point  The  global  model- 
in  principle-  distinguishes  on  a  finer  scale  the  aerosol  mixtures  over  various  regions-  The 
aerosol  extinction  at  1-315  /im  is  based  on  an  objectively  determined  surface  visibility  or 
meteorological  range  at  0,55  provided  by  the  GADS  model  data-  the  aerosol  modePs 
relative  extinction  (i.e.*  extinction  at  1.315  fim  relative  to  0,55  /im}-  and  the  particular 
aerosol  model's  vertical  profile  of  relative  extinction.  The  visibility  or  meteort)  logical  range 
(MR)  is  related  to  the  extinction  at  0,55  /rm  by  the  relation‘s 

MR  3.0/cr  (km)-  (4) 

where  o  is  the  surface  extinction  coefficient  at  0.55  //m.  Figure  3  shows  the  average  cam¬ 
paign  visibilities  for  several  locations  and  seasons  as  determined  from  GADS  climatology- 
Figure  3  illustrates  the  considerable  variation  between  regions  in  climatological  surface 
visibilities.  Here  Korea  exhibits  the  lowe-st  visibilities  (largest  visible  extinction)  while 
locations  in  the  United  States*  on  average-  experience  the  highest  visibilities- 
Figure  4  displays  a  typical  vertical  profile  of  total  absorption  and  extinction  based  on 
the  combination  of  HITRAN  and  the  GADS  model  using  input  from  a  thermosonde  mea¬ 
surement  at  Holloman  Air  Force  Base*  New  Mexico-  which  is  adjacent  to  the  White  Sands 
Missile  Range  (WSMR),  Note  that  the  surface  elevation  is  approximately  1,3  km  above 
MSL.  The  principal  variable  molecular  constituent  influencing  absorption  is  water  vapt>rnr 
humidity.  Near  the  surface-  the  molecular  line  absorption  at  1 .3 1 5  /r  m  is  typically  dominated 
by  water  vapt>r-  Every  significant  feature  can  be  associated  with  a  water  vapor  transition.  Up 
to  around  2  km  above  ground  level  (AGL)*  near  the  lop  of  the  atmospheric  boundary  layer, 
the  aerosol  scattering  magnitude  dominates  molecular  absorption.  Above  the  atmospheric 
boundary  layer,  line-by-line  absorption  is  the  dominant  contribution  to  both  absorption 
and  extinction.  Of  the  two  parameters-  molecular  and  aerosol  optical  depth,  the  aerosol 
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Fig.  4.  Vertical  profiles  of  total  absorption  and  extinction  at  1 .3 15  computed  using  the 
HITRAN  molecular  absorption  and  GADS  aerosol  databases.  These  profiles  were  computed 
using  data  from  a  thermosonde  profile  at  Holloman  Air  Force  Base.  New  Mexico,  on  6  June 
1 999  at  0208  MDT,  The  influence  of  the  boundary- layer  aerosol  in  the  extinction  curve  and, 
to  a  much  lesser  extent,  in  the  absorption  curve  is  evident  below  2  km  AGL. 

component  is  subject  to  more  uncertainty  due  to  its  reliance  on  climatological  means  and 
the  inability  to  directly  infer  its  vertical  distribution. 

Aerosol  absorption  (and  conduction)  contributions  are  also  dependent  on  the  HEL  in¬ 
tensity  magnitude^  ^  and  aerosol  vaporization.  As  discussed  in  detail  by  Sprangle  el 
the  aerosol  absorption  contribution  to  thermal  blooming  for  hygroscopic  aerosols  should  be 
normalized  by  the  term  (1  +  /?  ),  where  rjisa.  constant  of  order  unity  representing  the  ratio  of 
laser  energy  going  into  vaporization  to  laser  energy  conducted  into  the  air.  The  ratio  rj  is  not 
accounted  for  in  this  paper.  The  aerosol  absorption  coefficient  is  not  decreased  significantly 
by  the  presence  of  non  hygroscopic  aerosols.  For  a  moist  humid  lower  boundary  layer,  such 
as  found  near  the  Persian  Gulf  or  in  southern  New  Mexico  during  summer,  the  aerosol  ab¬ 
sorption  contribution  can  be  much  less  than  10%  of  total  absorption  at  1 .3 15  ^m.  Fora  dry 
boundary  layer,  such  as  in  Korea  in  winter,  the  aerosol  absorption  contribution  at  1 .3 15 
can  exceed  50%.  Above  the  atmospheric  boundary  layer,  the  aerosol  contribution  to  total 
absorption  is  typically  less  than  10%.  Prevailing  meteorological  conditions,  location,  and 
season  determine  the  relative  molecular  and  aerosol  absorption  contribution  as  a  function 
of  altitude. 


4.  Data  Sets  Analysis  with  Example 

Twenty-one  thermosonde  data  sets  from  different  geographical  locations  and  sea.sons 
were  used  in  the  analyses.  Their  location,  season,  and  number  of  usable  profiles  for  distortion 
number  calculations  are  listed  in  Table  2.  These  sets  or  “campaigns,”  each  consisting 
initially  of  approximately  20-30  vertical  profiles,  are  used  to  calculate  the  vertical  profiles 
of  molecular  and  aerosol  absorption  and  scattering  from  the  temperature,  pressure,  and 
humidity  profile  data,  as  well  as  to  determine  the  path-normal  transverse  flow  speed  using 
the  wind  velocity  measurements.  After  all  atmospheric  variables  in  a  vertical  profile  are 
screened  and  determined  to  be  suitable  for  analysis,  the  profiles  are  individually  analyzed 
for  up  to  84  slant  paths  <7  surface  distances  at  5-km  increments  for  12  upper-end  altitudes 
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Table  2.  Summary  of  range  and  mean  values  for  ainiosptieric  distortion  number  by 
thermo  sonde  campaign"" 


No,  of 

Year  files  Location 

Season 

N, 

(360  deg) 

Mean 

Nd 

(45  deg) 

Mean 

(90  deg) 

Mean 

1997 

16 

Bahrain 

Spring 

<20^200+ 

78,4 

<10-25+ 

21.9 

<10-20+ 

14.6 

1997 

26 

Bahrain 

Summer 

<20-150+ 

69.2 

<15-30+ 

22.3 

<10^20+ 

15.3 

1997 

18 

Riyadh,  SA  Fall 

<20-300+ 

92,1 

<10-40+ 

24.2 

<10-25+ 

14.9 

1999 

11 

Riyadh,  SA  Winter 

<10-150+ 

41.1 

<5-15+ 

12.5 

<5-15+ 

10,6 

1998 

6 

Riyadh,  SA  Spring 

<  1 0—  1 50+ 

57.6 

<10-20+ 

14.9 

<.5-15+ 

9.8 

1997 

15 

Korea 

Spring 

<20-200+ 

57.8 

<15-40+ 

24.9 

<10-35+ 

18,3 

1997 

11 

Korea 

Summer 

<30-200+ 

89.6 

<10-40+ 

28.1 

I9J 

1997 

14 

Korea 

Fall 

<20-200+ 

66.6 

<  10-40+ 

22.4 

<10-30+ 

15,3 

1998 

23 

Korea 

Winter 

<10-100+ 

37.7 

<5-20+ 

14.4 

<5-15+ 

9.8 

1998 

13 

Korea 

Spring 

<20-200+ 

87,7 

<10-40+ 

21.6 

<  10-30+ 

16.4 

1998 

19 

Australia 

Winter 

<10-200+ 

52,9 

<5-25+ 

18J 

<5-25+ 

15.1 

1998 

13 

Peru 

Winter 

<2C+200+ 

96.9 

<15-50+ 

31,6 

<10-40+ 

25,5 

1999 

19 

Korea 

Fall 

<1C+150+ 

54.6 

<10-25+ 

15.4 

<5-20+ 

11.1 

2000 

18 

Korea 

Winter 

<10-80-^ 

21.6 

<5-15+ 

<5-10+ 

5.7 

2000 

16 

Korea 

Spring 

<20-300+ 

79.9 

<  1 0-40+ 

22.9 

<10-30+ 

19,8 

1999 

19 

Qatar 

Fall 

<20-200+ 

56.2 

<10-25+ 

19.1 

<10-15+ 

13,8 

im) 

7 

Qatar 

Winter 

<10-150+ 

35,5 

<5-10+ 

8.1 

<5-10+ 

8, 1 

2000 

7 

Qatar 

Spring 

<10-150+ 

39.3 

<5-5+ 

6.8 

<5-5+ 

5,9 

1997 

10 

Holloman 

Fall 

<30-400+  194,8 

<15-50+ 

US 

<HF40+ 

19,5 

1998 

9 

Holloman 

Spring 

<10-200+ 

43,8 

<5-15+ 

12.0 

<5-15+ 

10.2 

1999 

20 

Holloman 

Summer 

<20-400+ 

85.9 

<10-35+ 

23.6 

<5-15+ 

18,9 

■*In  each  column,  ihe  largest  and  smallest  values  of  distortion  number  arc  highlighted  in  boldface.  SA. 
Saudi  Arabia. 


at  5()0‘m  vertical  increments).  Using  a  campaign  set,  the  means  were  computed  for  each 
individual  slant  path  for  each  campaign  and  plotted. 

The  most  frequent  reason  for  rejection  of  a  campaign  data  profile  was  missing  wind 
data.  For  “head-om*  (360  deg)  targets,  the  availability  of  usable  wind  data  was  especially 
irnptmam  for  slant  paths  that  extended  relatively  long  distances  at  low  altitudes.  Three  of 
the  campaigns  depicted  in  Fig,  3  (WSMRQO,  FIAFB_WAKE,  and  NH^WAKE)  were  not 
used  for  distortion  number  calculations  because  of  lack  of  good  wind  data  or  insufficient 
number  of  usable  wind  profiles.  In  all  3 10  ihermosonde  profiles  were  found  to  be  suitable 
to  be  analyzed  for  the  2 1  campaign  data  sets,  A  total  of  24,570  separate  slant  paths  were 
drawn  for  the  analysis. 

Figure  5  depicts  the  pattern  of  mean  distortion  number  for  the  Bahrain  summer  1997 
campaign  as  a  function  of  laser  altitude  abtTve  MSL  and  earth  surface  range  for  two  different 
aspect  angles,  360  and  45  deg,  for  a  ftK’used  laser  beam.  A  total  of  16  thermosonde  profiles 
were  analyzed  to  produce  these  plots.  In  the  upper  figure  (360-deg  azimuth),  the  surface 
target  lies  in  the  direction  of  the  aircraft  heading,  and  the  beam  transverse  Row  due  to  air 
speed  thu,s  is  minimal.  As  a  result,  Aj  values  exceeding  100,  a  critical  value  for  the  laser 
with  pha.se  compensation  according  to  Fig,  2,  appear  for  altitudes  below  1.5  km  MSL  and 
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Fig,  5*  Contour  analyses  of  mean  distortion  number  for  360-deg  (top)  and  45-deg 
(bottom)  azimuth  as  a  function  of  upper-end  altitude  (kilometers  MSL)  and  eanh  surface 
range  (kilometers)  forthe  Bahrain  summer  1997  measurement  campaign.  The  laser  aperture 
radius  is  50  cm,  and  the  laser  power  is  10^  W, 

earth  surface  distances  exceeding  7  km.  However,  when  the  target  lies  at  a  45-deg  azimuth 
angle  (lower  figure),  the  transverse  component  of  flow  is  sufficiently  large  so  that  the  largest 
values  of  are  less  than  35  and  confined  to  the  shorter  ranges.  The  shorter  ranges  in  this 
case  are  less  influenced  by  slew  velocity.  In  the  case  of  the  larger  azimuth  angles  >45  deg, 
uncompensated  blooming  effects  may  still  be  significant  but  phase  compensation  could  be 
used  to  minimize  atmospheric  effects.  The  mean  distortion  number  patterns  observed  in 
Fig.  5  generally  apply  for  the  other  campaigns  as  analyzed  according  to  azimuth  angle. 

5.  Discussion  of  Results 

Table  2  lists  by  thermosonde  mea.surement  campaign  the  range  of  distortion  numbers 
and  the  average  of  distortion  numbers  for  each  azimuth  angle  over  all  altitudes  and  ranges. 
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Fig,  6,  Plois  of  mean  distortion  number  by  season  for  different  measuremeni  locations 
(360-deg  a/jmuth). 
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Fig.  7*  Plots  of  mean  distortion  number  by  season  for  different  measurement  locations 
(45 -deg  azimuth). 


o 

□ 


+  Bahrain  1997 
A  Saudi  Arabia  1997-199S 
O  Korea  1997-1998 
□  Qatar  1999  -  2000 
•  Korea  1999  *2000 
JC  Holloman  1997  - 1999 


O 

& 


O 

I 


X 

□ 


— T" 
Fall 


Winter 


Spring 


Summer 


Season 


- 1 - 

- 1 - 

- F - - - 

- ^ - 

A 

+  Bahrain  1997 

rs 

A  Saudi  Arabia  1997-1998 

O 

w 

X 

0  Korea  1997- 1990 

f 

D  Qatar  1999-2000 

•  Korea  1999  -  2000 

+ 

0 

X  Holloman  1997  *  1999 

i 

n 

X 

A 

a 

□ 

• 

1 

Fall 

Winier 

spring 

1 

Summer 

Season 


in  each  column,  the  largest  and  smallest  values  of  distortion  number  are  highlighted  in 
boldface  and  boldface  italics,  respectively.  Figures  6-8  depict  plots  of  the  mean  distortion 
numbers  of  Northern  Hemisphere  campaigns  for  ail  altitudes  and  ranges  by  season  for  each 
azimuth  angle. 

For  azimuth  angles  of  360  deg.  the  typical  range  of  mean  Nj  is  from  less  than  10  to 
greater  than  200.  For  45  deg.  the  typical  range  of  mean  is  less  than  1 0  to  greater  than 
40.  For  90  deg,  the  range  of  N^j  is  from  less  than  5  or  1 0  to  greater  than  30,  Thus,  the 
slew  transverse  flow  component  plays  a  very  large  role  in  modulating  the  typical  range 
of  distortion  numbers  encountered.  It  is  also  evident  from  Table  2  that  for  the  Northern 
Hemisphere,  the  smallest  values  of  distortion  number  tend  to  occur  in  winter  and  the  largest 
values  tend  to  occur  in  the  spring  or  summer  seasons.  This  is  illustrated  in  Fig.  6,  which  i,s 
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Fig.  8.  Plots  of  mean  distortion  nymber  by  season  for  different  measurement  locations 
(90-deg  azimuth). 

a  plot  of  mean  distortion  number  by  season  and  geographic  location  for  the  360-deg  target 
azimuth*  Note  that  there  tends  to  be  a  greater  range  of  the  mean  values  of  distortion  number 
during  spring  than  in  other  seasons.  This  same  behavior  holds  also  for  the  azimuth  angles 
of  45  and  90  deg,  w^hich  are  depicted  in  Figs.  7  and  8. 

In  Eq*  (2),  the  distortion  number  is  dependent  on  the  absorption  coefficient  a.  Consid¬ 
ering  that  the  major  contribution  to  absorption  at  1315  /im  is  water  vapor  molecular  line 
absorption  for  altitudes  below  6  km.  a  measure  of  the  path-integrated  water  vapor  content 
.should  be  helpful  in  inferring  changes  in  the  distortion  number  between  seasons  and  locale. 
One  measure  is  slant  path  precipitable  water  (SPW).  This  represents  the  depth  of  water 
resulting  if  all  the  water  vapor  along  the  laser  path  between  the  two  ends  were  condensed. 
SPW  is  expressed  as 

SPW  =  f  dr  (cm),  (5) 

Jn  Pi 

where  w  is  the  water  vapor  mixing  ratio,  pj  is  the  density  of  dry  air.  and  p(  is  the  density 
of  liquid  water.  SPW  can  be  estimated  using  several  methods.  These  include  radiosonde 
measurements  of  humidity  (e.g..  Ref.  1 1 )  and  the  differential  absorption  lidar  (DIAL) 
technique  along  the  laser  propagation  path.'^  It  can  also  be  determined  using  passive  remote 
sensing  methtKls  such  as  microwave  radiometry^  or  refractive  path  delay  of  1 2- and  1 ,6-GHz 
carrier  signals  between  global  positioning  system  and  ground  stations.'® 

Figure  9  is  a  plot  of  mean  campaign  distortion  number  as  a  function  of  mean  SPW 
for  all  measurement  campaigns  listed  in  Table  2.  The  360-deg  azimuth  value  for  was 
used  in  the  comparison  since  this  angle  had  the  minimum  contribution  to  path  transverse 
flow  from  the  air-to-ground  slew.  The  mean  values  fora  campaign  include  all  air-to-ground 
paths  used  to  construct  Figs,  5-8*  generally  increases  with  increasing  SPW  up  to  about 
20  cm.  where  the  maximum  value  of  approximately  100  is  obtained  using  a  second- 
order,  least-squares  fit.  Leveling  off  is  indicated  when  SPW  exceeds  20  cm.  The  leveling 
may  be  a  result  of  increased  path  extinction  associated  with  the  increased  water  vapor  line 
absorption.  The  outlying  mean  value  Nfi^  194  is  associated  with  the  Holloman  fall  1997 
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SPW  (cm) 

Fig.  9.  Mean  atmospheric  distortion  number  vs.  mean  slant  path  precipitable  water  (centime¬ 
ters)  for  360-deg  a/Jmuthal  heading  for  campaigns  listed  in  Table  2.  Each  circle  represents 
a  campaign  mean  value.  Dashed  curve  marks  a  second-order,  least-squares  fit. 


campaign.  The  reason  for  such  high  magnitudes  is  probably  the  occurrence  of  very  low 
wind  speeds  during  this  particular  measurement  period  (Sept.  1997). 

6.  Summary  and  Conclusion 

Values  of  the  nondimensional  atmospheric  distortion  number  for  nonlinear  refraction 
have  been  calculated  for  different  locations  and  seasons  using  thermo  sonde  measurements. 
The  distonion  number  allows  inference  of  the  degree  of  nonlinear  refraction  or  thermal 
blooming  associated  with  an  atmospheric  HEL  path.  For  a  CW  laser  with  a  Gaussian- 
shaped  beam,  the  distortion  number  is  a  function  of  several  variables  including  the  HEL 
beam  characteristics,  atmospheric  variables  including  temperature,  moisture,  pressure,  and 
wind  velocity  and  laser  path  characteristics  including  location  or  position  and  slew  rate.  The 
laser  paths  used  in  the  calculations  encompass  a  range  of  air-to-ground  tactical  scenarios 
for  altitudes  up  to  6  km  ASL.  Scenario-dependent  calculations  are  accomplished  for  dif¬ 
ferent  geographical  regions  and  seasons  using  the  Air  Force  Research  Laboratory’s  global 
thermosonde  database,  the  HITRAN  molecular  absorption  line  parameters  database,  and 
OPAC/GADS  aerosol  extinction  models  for  the  COIL  wavelength  of  1.315  Subsonic 
beam  transverse  flow  is  assumed. 

Our  results,  using  fixed  HEL  parameters,  exhibit  a  wide  range  of  distortion  number  values 
depending  on  the  particular  path  scenario  and  azimuth  or  aspect  angle.  This  range  is  directly 
related  to  the  path’s  location  in  the  atmosphere  and  the  transverse  flow  associated  with  its 
slew  velocity.  For  the  case  of  head-on  azimuth,  distortion  number  values  can  exceed  100 
for  air-to-ground  paths  exceeding  10  km  in  length  and  originating  below  1 5  km  AGL.  The 
availability  of  accurate  wand  data  is  critical  to  this  case’s  determination.  For  paths  with  the 
larger  azimuth  angle  (>45  deg),  typical  values  did  not  greatly  exceed  40.  while  distortion 
values  at  the  largest  azimuth  angle  (90  deg)  were  even  lower.  Averages  of  the  distortion 
number  for  all  laser  scenario  paths  by  season  and  location  were  lowest  in  winter,  while  the 
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greatest  range  of  distortion  number  values  occurred  during  spring.  U  was  observed  that  the 
average  distortion  number  for  all  paths  during  a  measurement  campaign  correlated  with 
the  average  water  vapor  content  found  along  the  path.  This  correlation  is  primarily  due 
to  the  large  dependence  of  molecular  absorption  at  1315  fim  on  water  vapor  or  relative 
humidity. 

Considerable  uncertainty  exists  for  these  results;  thus  the  term  “estimate”  is  used  in 
the  title  of  this  paper  This  includes  the  inability  to  account  for  aerosol  absorption  and 
extinction  properties  as  a  function  of  the  HEL  irradiance  magnitude  by  neglecting  the  ratio 
of  aerosol  vaporization  energy  to  the  aerosol  energy  conducted  into  the  air,  the  uncertainty 
in  aerosol  extinction  associated  with  the  parameterization  of  aerosol  optical  properties  and 
structure  using  weather  and  aerosol  climatologies,  and  the  consideration  of  only  CW  and 
not  pulsed  laser  systems.  Also,  these  results  do  not  include  scintillation  eifects  due  to 
refractive  turbulence,  which  can  have  a  large  effect  on  the  laser  irradiance  field,  especially 
in  the  lower  atmosphere  regime  investigated  here.  The  results  may  prove  helpful  to  future 
research  efforts  that  seek  quantification  of  atmospheric  effects  on  HEL  propagation.  The 
same  analysis  approach,  if  proved  worthy,  could  be  applied  to  the  larger  set  of  global 
synoptic  rawinsonde  data  for  a  more  extensive  analysis. 
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